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O Abstract 

oo : 

We propose an extension of the cyclic hardening plasticity model formulated by 
-H- | Armstrong and Frederick which includes micropolar effects. Our micropolar exten- 

sion establishes coercivity of the model which is otherwise not present. We study 
then existence of solutions to the quasistatic, rate-independent Armstrong- Frederick 
model with Cosserat effects which is, however, still of non-monotone, non-associated 
type. In order to do this, we need to relax the pointwise definition of the flow 
rule into a suitable weak energy-type inequality. It is shown that the limit in the 
Yosida approximation process satisfies this new solution concept. The limit func- 
tions have a better regularity than previously known in the literature, where the 
original Armstrong-Frederick model has been studied. 

1 Introduction 

One of the well-known models to describe cyclic plasticity is the non-linear kinematic 
hardening model formulated by Armstrong and Frederick [3] . This model has been highly 
rated, because it is based on a physical mechanism of strain hardening and dynamic 
recovery, and because it has the capability of representing reasonably well the shapes 
of stress-strain hysteresis loops, especially those of constant strain ranges. Therefore, 
implementation of the Armstrong-Frederick model in finite element methods has been 
examined in several studies to date. Thus, that model is now available as a material 
model of cyclic hardening plasticity in commercial, general-purpose software for finite 
element analysis. 
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The Armstrong-Frederick model (AF) is a modification of the Melan-Prager model, 
which is well know in the literature and it can also be seen as an approximation of the 
Prandtl-Reuss model. The key modification of this simple model is the so-called " recall" - 
term, changing the evolution law for the symmetric backstress tensor b from a classical 
linear kinematic hardening law (Melan-Prager) to a nonlinear kinematic hardening law, 
i.e., 

b t = ce\ — d\e p t \b 

lin. kin. hardening recall-term, nonlinear hardening 

Here, e p is the symmetric plastic strain tensor, c and d are positive material constants. 
The "recall" -term entails the L°°-boundedness of the backstress b, a property which is an 
experimental fact since to the contrary, for high frequency cycles softening and rupture will 
occur. Therefore, the AF-model shows nonlinear kinematic hardening, but only to within 
a certain extent. The more realistic description of cyclic hardening plasticity experiments 
with the AF-model, however, has a prize to be paid: the model is non-coercive (bounded 
hardening), it is of non- monotone type and not of gradient type (non-associated flow rule). 
Thus, the AF-model is one of the prominent small strain plasticity models which has yet 
escaped the efforts of mathematicians to establish well-posedness. 

The mathematical analysis being quite challenging, there are no encompassing exis- 
tence results for this model in the literature. The first (partial) mathematical result for 
the Armstrong- Frederick model was obtained by the first author in the article [9]. There, 
the non-monotone, non-associated AF-model was written as a model of pre-monotone 
type (for the definition we refer to [2]). In this work the author used a Yosida approxima- 
tion to the monotone part of the flow rule. The obtained a priori estimates are, however, 
not sufficient to pass to the limit with such approximations and to get L 2 -strong solutions 
(see Section 3 in j9]). It was only shown that the limit functions satisfy the so called 
"reduced energy inequality". In the article [18] a regularization of the "recall" -term in 
the equation for the backstress was proposed. The existence of a rescaled in time solution 
to the Armstrong-Frederick model with the regularized equation for the backstress could 
then be established. The rescaling idea is very simple: a new time variable r = £(£) is 
proposed. Then the new system is easier to analyze, because the plastic strain is now 
uniformly Lipschitz with respect to the rescaled time. The main problem is to get back 
to the original system with the rescaled in time solution. It is, in principle, possible for 
rate-independent models, where the flow rule is invariant under scaling of the time vari- 
able. The Armstrong-Frederick model is rate-independent but the authors of [18] are not 
able to get back to the original system. The rescaling idea has already been applied in 
the plasticity context in [U [15], [16] . 

In this paper we want to extend the system of equations proposed by Armstrong and 
Frederick to include micropolar effects. In the classical metal perfect plasticity models 
at infinitesimal strains it has been shown in a series of papers [101 dH d2l E3] that a 
coupling with Cosserat elasticity may also regularize the ill-posedness of the Prandtl- 
Reuss plasticity model. This is possible because the Cosserat coupling leads to coercivity. 
Perfect plasticity, however, is yet characterized by a monotone flow rule of gradient type 
(associated plasticity). Therefore, the question arises naturally, whether adding Cosserat 
microrotations to the AF-model is still enough to regularize the problem in the way to 
satisfy the flow rule in a standard pointwise sense. From a modelling perspective, adding 
microrotations means to consider a material made up of individual particles which can 



rotate and interact with each other [201 EH (221 [25]. For phenomenological polycrystalline 
plasticity adding Cosserat effects is arguably a physically motivated regularization: the 
individual crystal grains are rotating and interacting with each other. 

The extension of the Armstrong-Frederick model to include Cosserat effects follows 
the lines proposed in [10], where the authors added the Cosserat effect to the classical 
elastoplasticity model with a monotone flow rule. There, and in our present approach, 
only the elasticity relation is augmented with Cosserat effects, the plastic constitutive 
equations, and notably the "recall" -term, is left unchanged, contrary to [15] . Regarding 
the effect of the Cosserat-modification for classical plasticity models, it has been proved 
that the new model is thermodynamically admissible and that there exists a unique, 
global in time solution to Cosserat elasto-plasticity. In [13] a H} oc - regularity result for 
the stresses and strains was proved, cf. [26j. The dynamic Cosserat plasticity was studied 
in [H], see also [121 E3]. Another way to regularize classical plasticity models is by 
introducing gradient plasticity effects [171 EH [271 GEJ EH] . However, a modification of the 
AF-model to include higher gradients will be left to future work. Moreover, the coupling 
with thermal effects can be treated as another attempt to regularize models from the 
inelastic deformation theory (cf. [HE]). 

Many non-monotone models from the theory of inelastic deformation processes in 
metals are also non-coercive (for the definition see [2]) and the existence results for such 
models is only very weakly examined. For example: the solutions obtained in articles [9], 
and [18] had a low regularity with respect to time and space (see also [301 EI], where the 
non-monotone model of poroplastisity was considered). In our opinion, it is expedient to 
first consider non-monotone but coercive models and to obtain better regularity results for 
the solutions. This article presents the first mathematical result in this respect for our new 
Armstrong-Frederick model with Cosserat effects, which is non-monotone, non-associated 
but coercive. 

2 The Armstrong- Frederick model with Cosserat ef- 
fect 

This section is devoted to the formulation of the Armstrong- Frederick model with Cosserat 
effects. 

From the mechanical results for Cosserat plasticity (see for example [TD], [H]) we 
conclude that we deal with the following initial-boundary value problem: we are looking 
for the displacement field u : Q x [0, T] — > M 3 , the microrotation matrix A : Q x [0, T] — > 
so (3) (so (3) is the set of skew-symmetric 3x3 matrices) and the vector of internal variables 
z = (e p ,b) : Q x [0,T] — > <S| ev x S^ ev (e p is the classical infinitesimal symmetric plastic 
strain tensor, b is the symmetric backstress tensor and the space 5| ev denotes the set of 



symmetric 3 x 3-matrices with vanishing trace) satisfying the following system of equations 

div x T = -f , 

T = 2fi(e(u)-e p ) + 2fi c (skew(V x u)-A) + Xti(e(u)-e p )H, 

—l c Aj-axl (A) = /i c axl (skew(V z M) — A) , 

ef G dI K{b) {T E y (2.1) 

T E = 2/x(e(u) - e p ) + Xix{e{u) - e p )TL, 

b t = ce p — d\e p \b, 

where e(u) = sym(V x u) denotes the symmetric part of the gradient of the displacement. 
The above equations are studied for x G Q C M. 3 and t G [0,T], where Q C M. 3 is a 
bounded domain with smooth boundary dQ and t denotes the time. 
The set of admissible elastic stresses K(b(x, £)) is defined in the form 
K{b) = {T E G S 3 : \dev(T E ) - b\ < cx y }, where dev (T E ) = T E - |tr(T B ) • 1, a y is a 
material parameter (the yield limit) and 11 denotes the identity matrix. The function 
Ik(p) is the indicator function of the admissible set K(b) and dlK(b) is the subgradient of 
the convex, proper, lower semicontinous function Ijctb)- The function / : Q x [0, T] — > 
M. 3 describes the density of the applied body forces, the parameters fi, A are positive 
Lame constants (the elastic constitutive equation can be generalized in the obvious way 
to anisotropic case), \i c > is the Cosserat couple modulus and l c > is a material 
parameter with dimensions [m 2 ], describing a length scale of the model due to the Cosserat 
effects, c, d > are material constants. The operator skew(T) = |(T — T T ) denotes 
the skew-symmetric part of a 3 x 3-tensor. The operator axl : so (3) — > M 3 establishes 
the identification of a skew-symmetric matrix with vectors in IR 3 . This means that if 
we take A G so(3), which is in the form A = ((0, a, /3), (— a, 0, 7), (— /3, — 7, 0)), then 
axl(A) = («,£, 7). 

Notice again that the expression \s p \b is a perturbation of Melan-Prager model - if d = 
then we obtain the Melan-Prager linear kinematic hardening model. 
The system (12. ip is considered with mixed boundary conditions for the displacement: 

u(x, t) = gD(x,t) for i6To and t > 0, 

T(x,t) ■ n(x) = g^f(x,t) for x G T N and t > 0, (2.2) 

where n(x) is the exterior unit normal vector to the boundary dQ at the point x, T E , 
Tn are open subsets of dQ satisfying dQ — Td U f^, T^ fl Tn = and "^(Td) > 0- 
%2 denotes the 2-dimensional Hausdorff measure. The system (12. ip is considered with 
Dirichlet boundary conditions for the microrotation: 

A{x,t) = A D (x,t) for xedQ and t > 0. (2.3) 

Finally, we consider the system (12. ip with the following initial conditions 

e p {x,0) = e p '°{x), b{x,0) = b°{x). (2.4) 

The free energy function associated with the system ( 12. ip is given by the formula 

pi/)(e,E p ,A,b) = n\\e(u) - e p \\ 2 + /j, c \\skew(V x u) - A\\ 2 

+ ^(tr(e(«) - e p )) 2 + 2/ c ||V,axl(A)f + 1||6|| 2 , (2.5) 



where p is the mass density which we assume to be constant in time and space. The total 
energy is of the form: 

£( y e,e p ,A,b)(t) = / pip(£(x,t),e p (x,t),A(x,t),b(x,t))dx. 
Jn 

From Section 2 of the article [9] we know that the inelastic constitutive equation occuring 
in ( 12. ip is of pre-monotone type. Moreover, if initial data (e p '°, b°) G iSf ev x iSf ev , then any 
solution (e p (t),b(t)) (if there exists) belongs to <S| ev x <S| ev , because the right hand side 
of ( 12. ip ^ is a subset of <S| ev and the following system 

— (tre p ) = with tre p (0) = 0, 
at 

— {trb) + d\e p t \trb = with tr6(0)=0, 

possesses the unique solution (tre p ,trb) = (0,0). Inspired by the work [9], we propose to 
rewrite the set of admissible stresses in the form 

K = {(dev (T E ),b) G S 3 dcv x S 3 dev : | dev (T E ) + cb\ < a y }, 

where the constant a y is the same as in the yield condition and the inelastic constitutive 
equation can now be written in the equivalent form 
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3 Main result 

Here we define a notion of the definition of the solution for the system (12. ip . Next we 
formulate the main result of this paper. First, let us start with the definition of L 2 -strong 
solution of the system (12. ip . 

Definition 3.1 Fix T > 0. We say that a vector (u, A, T, e p , b) G W 1,oo (0, T; H X {Q\ R 3 ) x 
H 2 (tt;so(3)) x L 2 (tt;R 9 ) x L 2 (fi;^ cv ) x L 2 (fi;5| ev )) is L 2 -strong solution of the system 

div x T = -f , 

T = 2^e( y u)~e p ) + 2^ c (skew(y x u)-A) + Xti( y e( y u)-E p )H, 

—l c A x axl (A) = /i c axl (skew( V x u) — A) , 

/ \ / dev (2»(s(u) - eP)) \ 

+ {d\e p \b) E dI ^{ -16 ) 

if 

1. \s p t \beL°°(0,T;L 2 (n,Sl v )), 

2. | dev (2fi(e(u(x,t)) - e p (x,t))) -b(x,t)\ < a y for almost all (x,t) G Q, x (0,T), 

3. the equations and inclusion above are satisfied for almost all (x,t) G Q x (0, T). 



d_fe p \.f 

dt 



Next, we are going to define a weaker notion of solutions to the system ( 12. ip . We give a 

motivation of this definition. It will be the energy inequality combined with special test 

functions. 

Let us consider another convex set (which will be used as set of test functions further on) 

KT = {(dev (T E ), -U) G S 3 dev x S 3 ev : | dev (T E ) - b\ + ±\b\ 2 < a y }, 

where the constant a y is the same as in the yield condition. It is not difficult to see that 
for all L 2 -strong solutions we get (for details we refer to figure 1 and 2 of [H]) 

^ (dev (T E ) - dev (f E j\ + b t ( - - b + - b) > for all (dev (f E ),b) G K* , (3.1) 

where T E G L°°(0, T; L 2 (Q, S 3 )) is any stress tensor. From the theory of elasticity, we know 
that there exists a positive definite operator C _1 : S 3 — > S 3 such that C _1 Te it = e t — e\. 
Integrating (13. ip over Q x (0,t) for t < T we obtain 

/ / C^ 1 T Ett (x,r)T E (x, r)dxdT -\ — / / b t (x,r)b(x,T)dxdT 
Jo Jq c Jo Jn 

< et(x,r)(T E (x,r) — T E (x,r))dxdT + / / C~ 1 T E:t (x,T)T E (x,T)dxdT 

Jo Jq Jo Jn 

+ - / / bt(x,T)b(x,T)dxdr = / / Vti t (x,r)(T(x,r) — T E (x : r))dxdT 
c Jo Jn Jo Jn 

— 2/i c / / ( skew(V x -u(x, r)) — A(x, r) )skew(V x u t (x, r))dxdr 

bt(x,r)b(x,r)dxdT . (3.2) 
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Let us assume that div T E G L°°(0,T; L 2 (Q,M. 3 )) (now the trace of T E in the normal 
direction to the boundary dfl is well defined) and T E (x, t) ■ n(x) = Qn(x, t) on Y N C dfl. 
Integrating by parts in the first term on the right hand side of (I3.2p . using equation f |2 . 1 1) i 
and the boundary data we have the following inequality 

- / C^ 1 T E (x, t)T E (x, t)dx + /i c / |skew(V x w(x, t)) — A(x, t)\ 2 dx -\ / \b(x, t)\ 2 dx 

2 Jn Jn 2c J n 

- / C' 1 T E (x,0)T E (x,0)dx + fi c / Iskew^^x.O))- A(x,0)\ 2 dx+— / \b(x, 0)\ 2 dx 
2 Jn Jn 2c J n 

/ u t (x,r)f(x,T)dxdr+ / / ut(x, r)divT E (x, r))dxdr (3.3) 

o Jn Jo Jn 

/ / 9D,t(x,r)(T(x,r) -f E (x,r)) -n(x)dSdr+ / C~ l T E j(x,T)f E (x,T)dxdr 
Jo Jv D Jo Jn 

- / / bt(x : r)b(x,T)dxdr — 2fi c / / (skewCV x u(x,r)) — A(x,T))A t (x,T)dxdT. 
c Jo Jn Jo Jn v J 



< 
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Using (12. ip ^ to the last term of the above inequality and integrating by parts in the last 
term on the right hand side of (13.31) we obtain 

- / C^ 1 T E (x,t)T E (x,t)dx + fi c / |skew(V x w(x, £)) — A(x, t)\ 2 dx 
2 Jn Jn 

+ 2/ c / \Vax\(A(x,t))\ 2 dx + — / \b(x, t)\ 2 dx < - / C' 1 T E (x,0)T E (x,0)dx 
Jn 2c Jq 2 Jq 

+ He |skew(V x u(a;,0)) - A(x,0)\ 2 dx + — / \b(x, 0)\ 2 dx + 2/ c / |Vaxl(A(x J 0))\ 2 dx 
Jn ^ c Jn Jn 

+ I u t (x,T)f(x,T)dxdr+ / / u t (x, t) div T E (x,r))dxdr (3.4) 



o Jn Jo Jn 

t r ft 



+ / / 9D,t( x ' T )(T(x,T) -f E (x,r)) ■ n(x)dSdr + / / C l T E ^(x,T)f E (x,r)dxdT 
Jo Jt d Jo Jn 

+ — bt(x,T)b(x,T)dxdr + Al c / / 'Vaxl(A(x,T)) -n 8od(An > t{x,T))dSdr, 

c Jo Jn Jo Jan 

where the boundary integrals are defined in the sense of duality between the space 
m (dQ; M 3 ) and H~* (dQ; M 3 ) (see [lj for details) . 

Let us assume that for all T > the given data F, go, gN, Ad have the regularity 

/ G H\0, T; L 2 (Q; M 3 )), g D G ff^O, T; H^ (T D ; R 3 )), (3.5) 

^eff^O.Tjff-^r^jR 3 )), AjjG #*((), T;#§(8fi;5o(3))). (3.6) 



Additionally let us assume that the initial data (e p '°, b°) G L 2 (Q;Sl cv ) x L 2 (fi;«S^ ev , 



|6°(x)| < - and | dev (T E (x)) - b°(x)\ < a y for almost all x G Q, (3.7) 

where the initial stress T E = 2//(e(w(0)) - e p '°) + Atr(e(«(0)) - e p '°)Il G L 2 (fi;5 3 ) is the 
unique solution of the following linear problem 

dw x T°(x) = -f(x,0), 
— l c A x axl (A(x, 0)) = fi c a.x\ (skew (V x u(x, 0)) — A(x,0)), 
u ( x > Q )\r D = 9d(x,0) T°(x) ■ n\ VN = g N (x, 0) A(x,0)\ m = A D (x,0), (3.8) 

with 
T°(x) = 2/x (e(w(x, 0))-£ p,0 (x))+2 / u c (skew(V x ?/(x, 0))-A(x, 0))+A tr(e(u(z, 0))-e p '°(x))ll. 



Definition 3.2 (solution concept — energy inequality) 

Fix T > 0. Suppose that the given data satisfy Ii3.5\) and A3.6\) . We say that a vector 
(u,T,A,sP,b) G L^Tj-ff^fyR 3 ) x L 2 (ft;S 3 ) x # 2 (ft;so(3)) x (L°°(fy<S| ev )) 2 ) sofces 
the problem l[2l\) -lj^4\) if 

(u t ,T t ,A t ,e p t ,b t ) G L^O,!;// 1 ^! 3 ) x L 2 (Q;R 9 ) x H 2 (Q;so(3)) x (L 2 (fi; S d 3 cv )) 2 ), 

£ae equations (j2.ip i and (j2.ip ^ are satisfied pointwise almost everywhere on Qx (0, T) and 
/or a// test functions (f E ,b) G L 2 (0,T; L 2 (fi; S 3 ) x L 2 (fi;S 3 ev )) snc/i tnat 

(dev(r £ ),6)G/C*, divf B GL 2 (0,T;L 2 (fi,M 3 )), T £ ■ n = g N on T N x (0, T), 

i/ie inequality 

— / C^ 1 T E (x,t)T E (x,t)dx + fi c / |skew(V x n(x,£)) — A(x, £)| 2 da; 
2 in in 

+ 2l c f \Vax\(A(x,t))\ 2 dx + ^- / |6(x,t)| 2 dx < i / C^T^T^a;)^ 
in 2c ifj 2 i n 

+ He \skew(V x u(x,0)) - A(x,0)\ 2 dx -\ / |6(x, 0)| 2 dx + 2l c / |Vaxl(A(x, 0))| 2 dx 

in 2c in in 

+ / / u t (x,r)f(x,T)dxdT+ / / u t (x, r)divr E (a;, r))dxdr (3.9) 

io in io in 



+ 



/ / g D ,t(x,T)(T(x } r) -T e (x,t)) ■ n(x)dSdr + / / C 1 T Eyt (x,T)f E (x,T)dxdT 
io ir D io in 

+ - / / &t(x, t)6(x, r)dxdr + Al c / / Vaxl(A(x,r)) • n axl(A D)t (x,r))dS'dT 
c io in io ian 

zs satisfied for all t G (0,T) ; wnere Tg G L 2 (Q;S 3 ) and (u(0),A(0)) G iJ 1 ^;^ 3 ) x 
id 2 (f2;so(3)) are unique solution of the problem Ii3.8\) . 

Theorem 3.3 (Main existence result) 

Let us assume that the given data and initial data satisfy the properties, which are specified 
in \3. 5\) - Ii3.8\) . Then there exists a global in time solution (in the sense of Definition 
\3.2\) of the system Ii2. 1\) with boundary conditions Ii2.2\) . Ii2.3\) and initial condition \2.J$ . 

Notice that the solution defined above has a quite nice regularity. We even get that 
e v t G L 2 (L 2 ), which yields \s^\b G L 2 (L 2 ). Unfortunately, this information is still not 
enough to obtain L 2 -strong solutions. This paper presents the first existence result for 
the Armstrong- Frederick model with Cosserat effects and the new solution concept. 
The proof of Theorem 13.31 is divided into two sections. First, we use the Yosida Approx- 
imation to the maximal monotone part of the inelastic constitutive equation. Next, we 
pass to the limit to obtain a solution in the sense of Definition 13.21 



4 Existence for the Yosida approximation 

We use the Yosida approximation for the monotone part of the flow rule from (12. ip in order 
to get a Lipschitz-nonlinearity only in equation (12. lp ^. Hence, we obtain the following 



system of equations 

div^T" = -/, 

T v = 2n(e{u v ) - e p ' u ) + 2/i c (skew(V,tO - A") + Xti(e(u v ) - e p ' u )l , 

-l c A x .axl (A") = fi c axl (skew(V x w y ) - A") , 

*T = ^{IdevCTD-yi-a^ ^™"^ (4.1) 

z/ | dev(Tg) — o"| 

7| = 2/x(e(u") - e"'") + Atr(e(w v ) - e p '")ll, 

&£ = cS P t ' V -d\6 P t ' U \b 1 '. 

The above equations are studied for x G Q C M 3 and t G (0, T). z/ > and {p}+ = 

max{0, p}, where p is a scalar function. 

The system (14. ip is considered with boundary conditions: 



u u (x,t) = g D (x,t) for x G T D and t > 0, 

T u (x,t) ■ n(x) = gjsr(x,t) for x G IV and £ > 0, 

A"(z,t) = A D (x,t) for iGffl and £>0 (4.2) 

and initial conditions 

e p ' u (x, 0) = e p '°(x), W{x, 0) = b°{x). (4.3) 



Theorem 4.1 Fix T > 0. Suppose that all hypotheses of Theorem \3.3\ are satisfied. Then 
for all v > t/tere exists a unique L 2 - strong solution (in the sense of Definition \3. 1\) 

(u\T\A\e v >\b v ) G W^^rjfT^njR 3 ) x L 2 (ft;M 9 ) x # 2 (ft;so(3)) x (L 2 (fi; S| ev )) 2 ) 

satisfying the system Iji4-1\ ) with boundary conditions (JJfy and initial conditions l[4-3\ ). 

Proof: The proof is very similar to that of a related result of the Armstrong- Frederick 
model without Cosserat effects: see Section 4 of [9 J . It uses similar techniques. We provide 
a sketch below. 

Step 1: In the beginning we formulate two lemmas. 

Lemma 4.2 Fix T > 0. Assume that (u u , T u , A u , e p ' u , b v ) is L 2 - strong solution of the 
problem ( [^. 1\ ) and \b°(x)\ < ^ for almost all x G fi. T/ien /or all v > 

|6"(x,t)|<^ for a. e. (x,i) G ft X (0,T). 

The Lemma [4.21 implies that if |6°(x)| < § then we can modify the equation (I4.ip fi in the 
following form 

^cer-dlefWO, ( 4 - 4 ) 

where 

n(6H \lA * l» > 5- 

II is a projection on the convex set {b G <S| ev : |fe| < §}. 

9 



Lemma 4.3 Fix T > 0. Assume that (u u ,T u ,A u ,e p,u ,b u ) satisfies the equation ^4-4\) ^ n 
the L 2 - strong sense on Q x (0,T) and \b°(x)\ < | for almost all x G Q. Then for all 
v>0 

U(b u (x,t)) = b u (x,t). 

For the proofs of Lemma 14.21 and T4.3I we refer to Section 4, Lemma 1 and 2 in [9]. From 
the above lemmas we conclude that if we find L 2 -strong solution for the system (14. ID with 
the equation (I4.4p instead of (I4.ip fi then this solution will be also solution of the original 
system (j4.ip . Using this information we define the function G v : S 3 x S^ cv — > <S| ev x <S| CV 
by the formula 

( I{|dev(T)-c6|-cr y } +1 £g=|{ [ , 

G,(T,6)={ (4.5) 

[ I {| dev (T) - cfc| - a y } + ^gl^, + f {| dev (T) - c &| - a y } + n(fe) . 

Let us consider the system ( 14. ip with the equation ( J4.4D instead of ( 14.1[) «. The flow rule 
for the modified system can by written in the form 

(er,K) = G u (-pV {eP , b) r(e u ,e^,A\bn), (4.6) 

where the free energy function i/j u is in the form 

pi) v {e v ,e p ' v ,A\b v ) = ti\\e(u") - £ p ' u \\ 2 + /j, c \\skew{V x u ,/ ) - A"f 

+ ^(tr(e(iO - en) 2 + 2/ c || V,axl(^)|| 2 + ^\\b»\\ 2 . 

Denote by S v (t) the energy associated with the system (14.1 p 

E v {u v e\e^\A\b v ){t) = f f nj; ,/ (u ,/ (x,t),s u (x,t),£ p ' ,/ (x,t), A u (x,t),b"(x,t))dx. (4.7) 



A fundamental tool in our proof is the following property of the energy function which 
results from our Cosserat modification: 

Theorem 4.4 (coerciveness of the energy) 

(a) (the case with zero boundary data) 

For all v > the energy function fl^.Tp is elastically coercive with respect to Vu. This 

means that 3C E >0,\/ue H%(£1), V A G H^(n), V e p G L 2 (tt), V b G L 2 (ft) 



(fi) + H^IL 2 ^) 



£»(u,e,e p ,A,b) > C E (\\u\\ 2 m{n) + \\A\\ 2 m 

(b) (the case with non — zero boundary data) 

Moreover, 3 C E > 0, V g D , A D G H^(dtt), 3 C D > 0, V e p G L 2 (fi) ; V 6 G L 2 (fi), 

V«G H 1 ^), V A G H 1 ^) with boundary conditions u\ gn = go and A\ an = Ad it holds 

that 



E"(u, e, e p , A, 6) + C D > C £ (llull^n) + P| 



m(n) + ll^llL 2 (n) 
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For the proof of Theorem 14.41 we refer to the Theorem 3.2 of the article |1Q] . Moreover, 
Lemma 3 of [9] provides the following properties of G u : 



• The function G u generates a bounded nonlinear operator from L 2 (Q;S 3 x S^ cv ) into 
L 2 (n;S*xSl v ). 

• For all (T^b 1 ), (T 2 ,b 2 ) G L 2 (fi;5 3 x <Sf ev ) the inequality 

(G U (T\ b 1 ) - G U (T 2 , b 2 ), (T 1 , b 1 ) - (T 2 , b 2 )) l2 > _£_Ji±l! [\\T l - T 2 \\\ 2 + \\b l - b 2 \\\ 2 ) 



holds. 

By the first statement above we can conclude that the system (14.11) with the flow rule 
(14. 5 p possesses the linear self-controlling property (for the definition we refer to jS]) and 
the second one means that this flow rule belongs to the class CAA - the class of Lipschitz 
perturbations of monotone vector fields. 

Step 2: The reasoning from above gives us that the system (14. ip with the flow rule 
(14. 5 p is coercive in the sense of Theorem 14.41 and it belongs to the class CAi with self- 
controlling property. The statement of the Theorem 14.11 may be proved using the general 
theory developed for the class CM. with self-controlling property (see for example [TJE]). 
This theory was also used to study the existence of solutions to a model of poroplasticity 
with Cosserat effect, which has a similar structure as the system (12. ip (details can be 
found in Section 4 of [E]). However, we would like very briefly to present the theory, 
which was used for the Armstrong-Frederick model without Cosserat effects (we refer to 
j9] for details). Let k > be a positive real number and C k be the following cut function 

Consider the following sequence of problems 

div x T k = -f, 

T k = 2fi(e(u k ) - e p ' k ) + 2/i c (skew(V x u k ) - A k ) + Xti(e(u k ) - e p ' k ) E , 

-l c A x axl (A k ) = ii c axl (skew(V x u fc ) - A k ) , 

ef = I { | d ev(T!)-^} + ^§§^, («) 

T k = 2fi{e{u k )-e p ' k ) + Xti{e{u k )-e p ' k )ll, 

b k = cC k (\s p ' k \)^--dC k (\ef\Mb k ). 

Here we drop the superscript v > and write (u k , A k , T k , e p,k , b k ) instead of 

(u u ' k , A u > k , T u ' k , e p,u,k , b v,k ). The system (14. 8 p is considered with boundary conditions (14. 2 p 

and initial condition (14.31) . 



Theorem 4.5 (global existence for Lipschitz nonlinearities) 



11 



Assume that for all T > the given data has the following regularity 

/GC 1 ^,!];^;! 3 )), g D G C\[0,T}; ^(r D ;M 3 )) , 

g N G ^([O^iy-^r^lR 3 )), A D g C\[0,T\; Hi (dQ; so(3))) , 

(^•°,6°) G L 2 (n;5L)xL 2 (f];5L). 

Then, for all k > i/ie approximate problem ( f^.ffi ) with boundary conditions \4-<fy an d 
initial condition ( f^.3| ) /ias a global in time, unique L 2 - strong solution (u k , A k , T k , e p,k , b k ) 
with regularity 

(u k ,A k ,T k ) G C l {[0,T];H l {Q;R 3 ) x # 2 (ft;so(3)) x L 2 (ft;M 9 )), 

^ G ^([O,!];! 2 ^^)), fe fc G C a ([0,r|;L 2 (n;^ w ). 

Notice that the function C k (\e p ' |) is an L°°-function therefore, the system (J4.8P includes 
global Lipschitz nonlinearities only and for a proof of this statement we refer to Theorem 
3.1 of the article DE. 



Step 3: Our goal is to pass to the limit with k — > oo. To do this we use the energy 
method and the properties of the function G u . 



Theorem 4.6 Let us suppose that all hypotheses of Theorem \4-l\ hold. Let us denote by 
£ k {t) the total energy associated with the system fl^.gp 

£ k (e k ,e p ' k ,A k ,b k )(t) = f L\e{u k )(x,t) - e p ' k (x,t)\ 2 + fi c \skew(V x u k (x,t)) - A k (x,t)\ 2 

+ ^(tr(e(u k {x,t)) - e M (x,t))) 2 + 2/ c |V x axl(A fc (£,£))| 2 

+ — \b k (x,t)\ 2 )dx. 
2c / 

Then for all t G (0, T) the following estimate 

£ k {e k ,e p ' k ,A k ,b k ){t) <C(T) 

holds and C(T) does not depend on k > (it depends only on the given data and the 
domain). 

Proof of Theorem \4-fy Calculate the time derivative of the energy and obtain 

j(s k (e k ,e p ' k , A k ,b k )(t)^j = J (2fi(e k - s p > k )(e k - e p ' k ) + Atr(e fc - e p ' k )ti(e k - e p ' k ) 

+ 2/i c (skew(Vw fc ) - A k ) (skew (\7u k ) - A k )\ dx + 4/ c f Vaxl(A*)Vaxl(A£) dx 

+ - f b k b k dx = [ T k Vu k dx- I dev (T k )ef dx - 2// c f (skew(Vw fc ) - A k )A k dx 
c Jn Jn Jn Jn 

+ Al c / Vax\(A k ) Vaxl(^) dx - / b k (ef - C k {\e p ' k \)^ T -)dx 
Jn Jn v \e P ' \s 

- f \dev(T k )-b k \\e p t ' k \dx-^ fc k (\ef\)\U(b k )\\b k \dx. (4.9) 



<o 
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Integrating by parts in the first and fourth term on the right hand side of (j4.9p . using the 
equations (14.81) i and (14. 8p ^ and boundary data we have (notice that ||^4|| 2 = 2||axM|| 2 ) 

£feV, «"*,*, *)(«)) < //«** + / m,,T*-n dS+ [ u* gN dS 
dt\ / Jn Jv D Jt n 

+ 4l c [ Vax\(A k )-na,x\(A D , t )dS- [ b k fef - C k (\e p ' k \)^-)dx. (4.10) 

J an Jn v \ £ t \ 

Integrating (14.101) with respect to time we have 

S k {e k , e p ' k , A k , b k ){t) < £ k {e k , e p ' k , A k , b k )(0) 

+ / [ f(r)u k (r)dxdr+ f I g D tt (r)T k : (r) • ndS 'dr 
Jo Jn Jo Jv D 

+ [ I u k (T)g N (r)dS + Al c f f Vaxl(A fc (r)) • n ax\(A Dtt (r))dS 
Jo Jv N Jo Jan 

~ f f b k {r)(ef{r)-C k {\ef{r)\)^^)dxdr. (4.11) 

Jo Jn v \St (r)! 7 

The continuity with respect to time yields that the initial values u k (0), A k (0) are solutions 
of the following linear elliptic boundary- value problem 

div x T fc (0) = -/(0), 

-l c A x axl{A k {0)) = fi c ax\{skew{V x u k {0))-A k (0)), 

u k (0) lrD =g D (0) T k (x) ■ n kN = g N (0) A k (0) lan = A D (0) , 

where 

T k (0) = 2fi(e(u k (0)) - e p '°) + 2^ c (skew(V ;c w fc (0)) - A k (0)) + Atr(e(w fe (0)) - e p '°)Il. 

This solution has the following regularity 

u k (0) e H'i^R 3 ), A k (0) e # 2 (ft;so(3)) 

and is independent on k. Hence, the initial energy S k (e k ,e p ' k ,A k ,b k )(0) is a constant. 
Integrating partially in time, we get 

f f{r)u k {T)dxdT = - I ! f t (r)u k (T)dxdr+ f f(t)u k (t)dx+ f f(0)u k (0)dx 
Jn Jo Jn Jn Jn 

< 2 J \\M T )Wh(nm dT+ 2j ll^^ll^^^^+ll/WII^^^H^WIU 2 ^;^) 

+ 11/(0)11^(^3)11^(0)11^(^3). (4.12) 

Applying Poincare's inequality to ||w ||i^ in ( 14. 12ft and the coerciveness of the energy 
function with respect to the gradient of the displacement vector we obtain 

/ [ f(r)u k (T)dxdr < C f £ k (e k ,e p ' k ,A k ,b k )(T)dT 

Jo Jn Jo 

+C\\f(t) \\ L 2 {nm (£ V, e p ' k , A k , 6 fc )(t)) 5 + C(t) , (4.13) 
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where the constants C, C(t) do not depend on k. In a similar way we show that 



o Jo, 



g N (r)u k t (T)dSdT < C £ k (s k ,e p > k ,A k ,b k )(T)dT 







+C\\g N (t) \\ H - i(rNm (£ k (e k , e p > k , A k , b k )(t)) § + C(t) . (4.14) 

The others boundary integrals in (14. lip are estimated using the trace theorems (the details 
will be shown in the proof of Theorem 15.11 of the next Section). It is not difficult to see 
that Lemma 14.21 holds for the sequence b k and 

/ yk( p,k r ik(\ p,k\\ E t \ j < I \uk\\ p,kij < II r fcn |lr p ' fc ll 

- / o [ft -<- {\ £ t \)TW\) ax ^ \ b W £ t \ax<\\b || L -(n;5| ev )lR IliH^iJ 

< C(H-a||^ fe ||| 2(n;5 3 ev) ) 

< ( equation (4.1) 4 ) < (4.15) 

< C(l + ^(112*11^^, + H&1i W a w) )) , 

where a is any positive number and C does not depend on k. By (14.1 2ft - ( I4TT51) we are 
able to obtain the following inequality 

£ k (e k ,e p ' k ,A k ,b k )(t) < a£ h (e k ,e p ' k ,A k ,b k )(t) 

+ C f £ k (e k ,e p ' k ,A k ,b k )(T)dr + C(T), (4.16) 

Jo 

where a is any positive number and the constant C, C(T) > does not depend on k > 0. 
Choosing a > sufficiently small and using Gronwall's lemma we complete the proof of 
Theorem 14.61 □ 

From the proof of the Theorem 14.61 we conclude that the sequence (u k , A k ,T k , e p,k , b k ) is 
bounded in W^°°{H l x H 2 x L 2 x L 2 x L 2 ). 

Step 4: To finish the proof of Theorem 14. II we need L°°(L 2 ) -strong convergence of the 
sequence (T fc , b k ). 

Theorem 4.7 Let us assume that the given data satisfy all requirements of Theorem \4-6] 
Then, 

£ k (e m - e n , e p ' m - e p ' n , A m - A n , b m - b n )(t) — ► 

for m,n —¥ oo uniformly on bounded time intervals. 

For the proof of the above Theorem we refer to Section 4 of [9J. Theorem 14.71 implies 
the L°°(I/ 2 )-strong convergence of the sequence (T k ,b k ) hence we can pass to the limit 
in the system (14. 8 ft with k — > oo and obtain L 2 -strong solution of the system (14. ip . The 
uniqueness follows immediately from coerciveness of the energy function evaluated on the 
difference of two solutions of the system (14. ip . The last statement finishes the proof of 
Theorem gj] □ 
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5 Proof of the Main Theorem 13.3 



In this section we are going to prove the L 2 (L 2 )- boundedness for the time derivatives of 
the sequence (u u , T u , A v , e p ' u , b u ). It is the main part to prove the Theorem 13.31 

Theorem 5.1 (Energy estimate) 

Assume that the given data and initial data satisfies $3.5\) - ( fff.ffj) . Then for all t G (0, T) 
the following estimate 



1 



{| dev {7%){t) - b»(t)\ - a y }\dx + / / CT^^T^r^r 



n 2^ Jo Jn 

-t r ft 



+ 2fic f f |skew(VX(^)) -A^{r)\ 2 dxdr + U c f f \Vajd{A u t {T))\ 2 dxdr < C{T) 
Jo Jn Jo Jn 

holds and C(T) does not depend on v > (it depends only on the given data and the 
domain). 

Proof: Compute the time derivative 

Jt(i^ {l deV ^W " ^*)l " a ^ dx ) = j/t U {t)(&ev{Tl t ){t) - W t (t))dx 

= (by the elastic constitutive relation) = / e v t {t)T Et {t)dx — / C _1 T^ t T Et {t)dx 

Jn Jn 



er(W(t)dx. (5.1) 

From Lemma 14.21 we have the following inequality 

/ e^ ,v (t)b v t {t)dx = (by the equation for the backstress) = (5.2) 

Jn 

c [ \e p /{t)\ 2 dx -d [ \e p /{t)\e?/{t)b v {t)dx > [ \e p t ' u (t)\ 2 (c - d \b u (t)\)dx > . 
Jn Jn Jn 

Notice that 

e u t {t)T v Et {t)dx = J Vu u t {t)T»{t)dx - f 2/i c (skew(VX(£)) - A$(tj\ skew {V x u v t (t))dx 

Vu v t {t)T»{t)dx - 2/i c f |skew(V a ,.<(t)) - A v t {t)\ 2 dx 
n Jn 

- 2/i c / (skew{V x u u t {t))-A u t {t))A u t {t)dx. (5.3) 

Using (I4.ip ^. integrating by parts in (15.31) we have 



e u t (t)T^ t (t)dx= / u v t (t)f t (t)dx + / g Dtt (t)(T t u (t))-ndS+ / u v t {t)g N , t {t)dS 

n Jn Jv D Jv N 

- 2fx c I |skew(V x <(t)) - A v t {t)\ 2 dx - 4/ c f \V&x\(A u t (t))\ 2 dx 
Jn Jn 

+ Al c I Vs^\{A u t {t))-ns^\{A D . t {t))dS. (5.4) 

Jan 
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The first term on the right hand side of (15.41) is estimated as follows 

I u v t {t)ft{t)dx < ||<(t)|| H i ( n^3)||/t(t)||^( ni R3) • (5.5) 

Jn 

We estimate the appearing boundary integrals 

/ 9nAt)mt)) ■ ndS < \\9nAt)\\ HhrDm \\Tm ■ n\\^ 

< (from the trace theorem in the space Lj iv (Q) see for example in [18]) < 

< C(||^(t)|| ia( n i5 8 ) + \\dwT t v (t)\\ LHnm )\\g D At)\\ Hi{m]R3) 

< C||^(t)||i^ (n ^-)I|to,*(t)|| j4(0OiRa) + C||/ t (t)I| £a( o iH a ) ||^ tt (*)||^ 4(ao ^, ) , (5.6) 

where C > does not depend on v. Moreover 

u v t (t)g Nt (t)dS < \\QNt{t)\\ i, ,JK(t)|| i 

J- JV 

< lbiv 1 t(*)|| i? _i (rjv . ffi 3 ) |K(*)||^i(n;M3). (5.7) 

We use the ^-regularity of the microrotations to estimate the following integral 

(Vaxl^i))) -naxl(A D At)) dS 

an 

< ||Vaxl(^(t))) • n\\ H ^ {dnm \\ a AA D Am Hhdnm 



< C ||Vaxl(^(t))|| L2(n . R9) + ||Aaxl(^(t))|| L2(n;R 3 ) ||axl(A 



w,t 



H2(an;i 



C||Vaxl(^(t))|| L 2 (n . R9) ||axl(A At (t))|| 



H?(an-; 



+ C^||skew(V<(t)) - A^t)\\ LHnm \\^\(A D At))\\ Hh{m;R3) • (5-8) 

Inserting (I5.2p - (l5.8p into (15. ip . using Cauchy's inequality with a small weight at the ap- 
proximate sequence and integrating with respect to time we obtain the following inequality 

/ i{| dev (T»)(t) - r(t)\ - a y }\dx + f [ C-^Tffl^dxdT 

Jn Av Jo Jn 

+ 2fi c f f\^ew(V x u u t (T))-A v t {T)\ 2 dxdT + Al c f f \Vaxl(A u t (r))\ 2 dxdT 
Jo Jn Jo Jn 

< aj^ (||ir(r)||i a(n ^ ) + ||skew(V<(r)) -^(r)||i 2(n;R9) + ||Vaxl(^(r))||i 2(f , ;M9) )dr 

+ a f K(r)\\ 2 H i ( a- m dr + C(T,a)+ f l{|dev(T£)(0) - 6"(0)| - a y }%dx, (5.9) 
Jo Jn Av 

where a > is any positive number and C(T) > does not depend on v (it depends 
on the given data and the time interval only). We observe that T u (0) e L 2 (Q;S 3 ) is the 
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unique solution of the problem 

div x T"(x,0) = -/(x,0), 
-l c A x axl(A u {x,0)) = n c 8ixl(skew(V x u u (x,0))-A v (x,0)), 
u v {x, 0)\ Fd = g D (x, 0) , V(x, 0) -n lrN = g N (x, 0), A x > Q (x, 0), an = A D (x, 0) , 

where 

T u (x,0) = 2^(e(u u (x,0))-e p '°(x)) + 2ij c (skew(V x u u (x,0))-A u (x,0)) 

which implies that dev (TJJ(O)) = dev (Tg) and b u (0) = b°. From the assumption (13.71) we 
get that the last term on the right hand side of (I5.9P is equal to zero. 
We know that the function ||Vw||.L 2 (n;]R 9 ) + f r \u\dS is a norm on i/^f^R 3 ) equivalent 
to the standard norm (see [32]). Applying this fact in (15.91) and the following well-known 
estimate [[IS], p. 36] 

l|Vw||£ 2(n) < C^ v rl (||div u\\ 2 L2{n) + ||curl u\\ 2 L2{n) ), 

(the constant C^J 1 does not depend on u and 'curl' is the rotation operator) we get the 
following inequality 

/ i{| dev (T»)(t) - r(t)\ - a y }\dx + f I C^T^T^dxdr 
Jn Av Jo Jn 

t r rt 

i/C \|2 



+ 2fi c \skew(V x u u t (T))-A u t (r)\ 2 dxdT + Al c \Vax\(A u t (r))\ 2 dxdT 

Jo Jn Jo Jn 

< a jT (\\T?(t)\\ 2 l2{u . s3) + ||skew(V<(r)) - 4"tollW) + l|VaxlK(r))||| W) )dT 

+ a J (||div<(r)||| 2(n . R) H-||curl<(r)||i 2(n . R3) )drH-C'(r). (5.10) 

From the observation div u v t = tr s(u") — tr e P,u we have 

f ||div <(r)||| 2(n . R) dr <C f \\T?(T)\\l Hn . s3) dT , (5.11) 

Jo Jo 

where the constant C does not depends on v. Now we choose a > sufficiently small and 
we arrive at the inequality 



1 



{| dev (T%)(t) - P(t)\ - a y } 2 + dx + / / C^Tl^T^dxdr 



n 2^ jo Jn 

+ 2/i e [ [\s\Lew{V x u v t {r))-A u t {r)\ 2 dxdT + Al c f f \V txx\{A v t {r))\ 2 dxdr 
Jo Jn Jo Jn 

< a I \\cw\u" t (T)\\ 2 L2{a . m dT + C(T,a). (5.12) 

Jo 
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Notice that 



2fi c J J |skew(VX(r)) - A u t (r)\ 2 dxdr + Al c f f \Va,xl(A^(r))\ 2 dxdr 
Jo Jn Jo Jn 

= 2^ c / \skew{V x u u t {r))\ 2 dxdT + 2fi c \A u t (T)\ 2 dxdT 

Jo Jn Jo Jn 

- AfjL e f /" skew(VX0")) • AZ{r)dxdr + 4l e [ [ |Vaxl(^(r))| 2 dxdr 
Jo Jn Jo Jn 

> \i c / / |skew(V x .<(r)) \ 2 dxdr-2/2 c / / \A v t {r)\ 2 dxdT 

Jo Jn Jo Jn 

+ 4/ c / / \Vaxl(A^(r))\ 2 dxdT > ( Poincare's inequality ) > 

Jo Jn 

> He |curl u u t (r))\ 2 dxdr - Afi c C n / / |V axL4£ (r)| 2 dxdr 

Jo Jn Jo Jn 



t 



+ 4/ c / / \Vax\(A»(T))\ 2 dxdT, (5.13) 

Jo Jn 

where the constant Cn > depends on the domain Q only. Shifting the negative term on 
the right hand side of (15.131) into the left we may obtain the inequality 

2/i c I [\skew(V x u»(T))-A»(T)\ 2 dxdT + 4l c [ f \Vax\(A u t (r))\ 2 dxdr 
Jo Jn Jo Jn 

> C(n,fj, c ,l c )(fj, c I I \cm\u v t (T))\ 2 dxdT + Al c I f \\/ajd(A u t (r))\ 2 dxdT\ (5.14) 
^ Jo Jn Jo Jn ' 

If we apply inequality (15.141) to (I5.12p then we can again choose a > sufficiently small 
to end the proof. □ 



Theorem 5.2 Let us suppose that all hypotheses of Theorem \5. 1\ hold. Then the sequence 
{(er,6?)}„>o is bounded m L 2 (0,T;L 2 (n;Sl v ) x L 2 (fi; S d 3 cv )). 

Proof: The proof of the Theorem 15.11 yields that the sequence {u^} u> q is bounded in 
L 2 (0, T; H l (Q; R 3 )). From the elastic constitutive equation we obtain 

l|eT(*) 11^(0^) < IK<(*))II*W) + \\C- l T»(t)\\ LHQ]S3) 

< \\e(u»(t))\\ LH n.,ss) + C\\1?(t)\\ v >(ci i s>), (5.15) 

where C > does not depend on v > 0. Theorem 15.11 gives us the boundedness of the 
sequence {st' U }v>o i n L 2 (0, T; L 2 (Q; «S| ev )). Next using equation (j4.ip fi we have 

VU\\\2 S h( £\\<-V> V ( + \\\1 I rfl II \r-P> V ( + \\hV( + \ l|2 



II^WIIi 2 (n ;5 3 cv ) < ^Kll^(*)lli» (n ^ r ) + d a ll|er(*)r(*)lli»(n ! 53 w) ) (5-16) 



and 



\cr(tW(t)\\l H n 



5 dov) 



er(t)\ 2 \b»(t)\ 2 dx 



< 



_p,u 



L 2 (.n-,s'i cv )\ 



\b u (t)\ 2 \\ L °°(n;Sl v ) 



(5.17) 



From Lemma [4. 21 we know that \b u (x,t)\ < % for almost all (x,t) G ft x (0,T), then the 
term on the right hand side of (15.171) is bounded for almost all t G (0,T). This finishes 
the proof. □ 



Proof of Theorem GO 

The motivation of the Definition 13.21 yields that L 2 -strong solutions of the system (14.11) 
with boundary and initial conditions (14. 2 p and (14. 3 p satisfy the energy inequality, which 
means that for all test functions {T E ,b) G L 2 (0, T; L 2 (ft; <S 3 ) x I? (ft; <Sf cv )) such that 

(dev(f B ),6) G/C*, div Tg; G L 2 (0,T; L 2 (ft,M 3 )) , f E ■ n = g N on T N x (0,T) 
the inequality 



1 



i-i™ 



(x)dx 



T E (x,t)T^(x,t)dx + fi c / |skew(V x ^(x,t)) -^(x,t)| 2 c/x 
+ 2/ c / |Vaxl(^(x,£))| 2 dx + ^- / \b v {x,t)\ 2 dx < J / d- 1 !^)^ 
+ He \skew(V x u(x,0)) - A(x,0)\ 2 dx -\ / \b(x, 0)\ 2 dx + 2l c / |Vaxl(A(x, 0))| 2 dx 

+ / / u^(x,r)f(x,T)dxdr+ / / m^ (x, r)divT,E(;r, r))dxdr (5.18) 

io in io in 

/ 9D,t{x,T){T v {x,T) —Te(x,t)) ■ n{x)dSdr + / / C^ 1 T Et (x,r)TE(x,r)dxdT 
o Jt d io in 

+ - bt(x,r)b{x,T)dxdT + 4l c / Vaxl(A 1/ (x,r)) • n asd(A D)t {x,r))dSdT 

c io in io ian 

holds for all t G (0, T), where Tg G L 2 (ft; <S 3 ) and (w(0), A(0)) G if 1 (ft; M 3 ) x # 2 (ft;so(3)) 
are unique solution of the problem ( I3.8p . From Theorem I5.1l and l5.2l we conclude that for 
a subsequence (again denoted by v) we have 

u u ^u in L°°(0,T;H\n;R 3 )), 

u u t ^u t in ^(O,!;^ 1 ^;! 3 )), 

A u -» A in L°°(0,T;# 2 (ft;so(3))), 

T V ^T in L°°(0,T;L 2 (ft;M 9 )), 

T^^T t in L 2 (0,T;L 2 (ft;M 9 )), 

F-& in L~(0,T;L 2 (ft;<S 3 ev )), 

W t ^b t in L 2 (0,T;L 2 (ft;S 3 cv ))- 
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The right hand side of (I5.18P has the limit for v — > + , the left one is a convex functional 
thus, the sequential weak lower semicontinuity implies that 

- / C^ 1 T E (x, t)T E (x, t)dx + /i c / |skew(V x w(x, t)) — A(x, t)\ 2 dx 
2 Jn Jn 

+ 2/ c / \Vax\(A(x,t))\ 2 dx + — / \b(x,t)\ 2 dx 
Jn 2c J n 

< liminfi- / C~ l T E (x, t)T E (x, t)dx + \x c I \skew (V x u u (x, t)) - A u (x,t)\ 2 dx 

u^q+ y 2 j Q - - j Q 

+ 2l c [ \Vax\(A u (x } t))\ 2 dx+^- [ \b u (x,t)\ 2 dx\. (5.19) 

Jn 2c J n J 

Using (15.191) we complete the proof of the Theorem 13.31 □ 
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